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1 •  Introduction 

A  3  dimensional  (i.e.,  3-factor)  design  involves  3  factors 
treatment,  row  and  column.  The  sets  of  treatments,  rows  and  columns 
are  the  sets  of  levels  of  3  factors.  The  determination  of 
connectedness  of  3  and  higher  dimensional  designs  is  a  difficult 
task.  The  method  of  Srivastava  and  Anderson  (1970)  involves  the  use 
of  X-operator.  The  method  of  Birkes,  Dodge  and  Seely  (1976)  known  as 
BDS  algorithm,  uses  the  incidence  matrix  of  a  design  to  find  the 
estimable  parametric  functions  by  R-and  Q-processes.  The  X-operator 
is  indeed  a  powerful  mathematical  tool  and  R-  and  Q-processes  are 
indeed  good  mathematical  algorithms.  When  it  comes  to  actual  doing 
and  finding,  we  find  those  two  earlier  known  methods  are  difficult  to 
implement  and  there  is  something  conspicuously  missing.  Our  search 
leads  to  a  new  concept  of  connected  path  for  three  dimensional 
designs,  a  generalization  of  a  powerful  idea  of  Bose  (1947)  for  two 
dimensional  designs.  The  proposed  method  takes  into  account  the  idea 
of  reduction  in  dimensionality,  discussed  in  Section  2,  and  the 
concept  of  connected  path,  discussed  in  Section  3,  and  some  graph 
theoretic  tools.  The  method  consists  of  drawing  graphs,  i.e.,  joining 
points  by  lines  directly  from  the  design  and  does  not  require  to  write 
down  the  incidence  matrix  and  construct  another  matrix  from  it  to 
perform  R-  and  Q-algorithmic  processes  like  in  BDS  algorithm. 
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2.  Reduction  in  Dimensionality 

Consider  a  3  dimensional  design  with  a  rows,  b  columns  and  v 
treatments.  The  standard  additive  model  for  the  observations  is 


E(yijk)  “  M  +  °i  +  6j  +  V  i"1»****a*  j“l, •••.!>,  k-l,...,v,  (1) 

where  the  observations  y,  's  are  correlated  or  uncorrelated  random 

ijk 

variables  with  equal  or  unequal  variances.  Consider  a  pair  of  columns 
the  row  i  and  the  pair  of  treatments  (k,k')  occurring  in  the 
design  corresponding  to  i  and  (1 ,j  ' ) .  Define  x  *  x  -  x  , 

j  K  K 

0  *»  0  -  0.  Then  we  have  from  (1), 

\J  *J  )  j  J 

E(yij’k'  "  yijk}  “  B(j,j  ’)  +  T(k,k’)’  J*j’m  (2) 

Similarly  for  a  pair  of  rows  (i,i'),  the  column  j  and  the  pair  of 

treatments  (k,k')f  we  have 

E(yi'jk'  "  yijk)  *  °(i ,i ' )  +  T(k,k' ) ’  1*i'’  (3) 

where  a.  , .  ■  a  ,  -  a  .  The  pair  of  columns  (or.  the  pair  of 

V,  1 » 1  ^  1  1 

rows  (i,i'))  can  be  treated  as  a  one-dimensional  block.  We  thus 


reduce  the  dimension  of  a  3  dimensional  design  to  two  by  considering 
all  pairs  of  columns  (or  rows)  and  all  pairs  of  treatments.  The 
equations  (2)  and  (3)  are  essentially  less  than  full  rank 
reparametrization  of  (1). 

Denote  the  sets  whose  members  are  the  pairs  of  rows,  columns  and 
treatments,  respectively  by 

I  -  {(!,!') ,  i*i\  i,  i*  e(l,...,a)}, 

j  -  j*j\  j,  .r  e(i . w}, 

K  ■  {(k,k'),  k,k'  e( 1 , . . . ,v)  } . 


(4) 


The  cardinalities  of  the  three  sets  are  jlj  =  a(a-l),  |jj  =  b(b-l)  and 
| k |  “v2.  In  case,  the  treatments  are  not  replicated  more  than  once 
in  a  row  or  a  column  of  the  design,  we  assume  k*k'  in  K.  Then  |k|  = 
v(v-l)  . 

We  write  k~k*  if  two  treatments  k  and  k*  are  connected,  i.e., 

t  is  estimable.  Similarly  for  rows  and  columns. 

(k,k ' ) 

Definition  1.  Two  members  (k^  ,k2)  and  (k3  ,k4)  in  K  are  said  to  be 
connected,  (kltk2)  -  (k3,k4),  if  T(k3>k4)  “  T(k1,k2)  ls  estimable 
under  (1). 

Note  in  Definition  1,  the  cases  k1=k2,  k3=k4  and  kj=*k3  ,  k2=k4  are 

trivially  true.  Clearly,  ~  is  an  equivalence  relation  on  the  set  K. 

Let  K  ,  K  , . • . , K  ,q(>  1)  is  an  integer,  be  the  distinct  equivalence 
12  q  — 

classes  of  ~  on  K. 

Lemma  1.  If  (k,.kj  ~  (k„,k,)  then  (k  ,k  )  ~  (k.  ,k.  ),  where 
-  1  Z  j  4  uj  U2  U3  U14 

(uj  ,u2  ,u3 ,u4)  is  any  permutation  of  (1,2, 3, 4)  in  the  dihedral  group  D4 
of  permutations  of  order  8. 

Lemma  2.  The  following  are  true. 


a. 


(kj^)  ~  (k2>k3)~. ..^.(k^  ^k^) 

where  w(>^  2)  is  an  integer, 


(k  ,k  )  =->  k  -~k_k  _ _ „k  ~k  , 

w’  1  12  3  w-1  w’ 


b.  kp-k2  •“>  (kj  ,k)  ~  (k2  ,k)  for  any  k, 

c.  (kj ,k2 )  ~  (k3,k4)  and  k3  ~  k4  — >  kj  ~  k2  , 

d.  kpk2 ,  kj~k4  (k^,k2)  ~  (k3,k4). 

The  proofs  of  Lemmas  1  and  2  are  easy.  Lemmas  1  and  2  are  very  useful 

in  finding  the  equivalence  classes.  To  illustrate  this,  suppose 
(kj  ,k2 )  ~(k3,k4),  (k2,k5)  ~  (k4,k6)  and  (k5,k6)  ~(k6,k5).  Then, 


•  r 


by  Lenina  1,  we  have  (ki,k3)  —  (k2,k4)  and  (k2,k4)  ~  (ks.kg).  Thus, 

I 

(k1,k3)  ~  (k5,k6).  It  now  follows  from  Lemma  2  that  kj~k3  and  k2~k4 . 

Note  that  ciT.,  ,,  .  .  ,  ci  and  c2  are  real  constants.  Is  a 

constrasts  of  t’s.  However,  In  Definition  1,  we  treat  (kj,k2)  and 

(k3,k4)  as  two  treatments  in  the  reduced  2  dimensional  design  and 

consider  a  simple  contrast  t„  ,  .  -  t,.  ,  v  (i.e.,  c,  •»  c9).  The 

(k3,k4)  (k i  ,k2)  1  i 

estlmability  or  nonestimabillty  of  ciT.,  ,  -  c9t.,  ,  .  can  easily 

1  (k3,k4)  <■  (k !  ,k2) 

be  determined  from  the  equivalence  classes  ,K^ , • • • ,K^ . 

Theorem  1.  A  3  dimensional  design  is  connected  w.r.t.  treatment  if 
and  only  if  (iff)  any  two  members  in  K  are  connected. 

Proof .  Suppose  any  two  members  in  K  are  connected.  Then,  for  an  two 
treatments  kj  and  k2,  we  have  (kj  ,k2)  ~  (k2,k1).  Thus  kj~k2  by  Lemma 
2.  If  the  design  is  connected,  then  for  any  two  pairs  (k3,k2)  and 
(k3,k4)  in  K  we  have  kpk2  and  k3~k4.  Therefore,  (k1,k2)  —  (k3,k4)  by 
Lemma  2.  This  completes  the  proof. 

Theorem  2.  A  3  dimensional  design  is  completely  connected  (i.e., 
connected  w.r.t.  treatment,  row  and  column)  iff  any  two  members  in 
each  of  I,  J,  and  K  are  connected. 

Proof .  Similar  to  the  proof  of  Theorem  1. 

If  the  design  is  connected  w.r.t.  treatment  then  q  -  1,  i.e.,  there  is 
one  and  only  one  equivalence  class.  In  case  the  design  is  not 
connected  w.r.t.  treatment  we  have  q  !>  2  and  any  two  members  in  an 
equivalence  class  are  connected  and  any  member  in  an  equivalence  class 
Is  not  connected  to  any  member  in  a  different  equivalence  class.  The 


.v.v.%  .•  .v.v 
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equi valence  classes  K  ,K  give  all  the  estimable  and  non- 

1  2  q 

estimable  row  and  column  contrasts  from  the  equivalence  classes 

I,,...,  I  and  ,  p,r(>  1)  Integers,  of  ~  on  I  and  J. 

1  r  1  p  — 

It  Is  easy  to  see  that  any  two  members  in  K  occurlng  In  two 
distinct  rows  (or,  columns)  of  (or,  (1,1'))  are  connectd.  Two 

pairs  of  columns  ( j  i » J 2)  an<*  ( J  3  »J  i+)  are  connected  if  a  pair  of  treat 
ments  in  ( J 1 » j 2 )  *s  connected  to  a  pair  of  treatments  In  ( j  3 , j  4 ) •  If 
(jl,j2)  and  ( j 3  > j 4 )  have  a  common  pair  of  treatments  then  (j1,j2)  ~ 

( j 3  » j 4  )  or,  in  other  words,  every  pair  of  treatments  in  ( J  1 » j 2 )  Is 
connected  to  any  pair  of  treatments  in  ( J  3 , j  4 )  •  Similarly,  two  pairs 
of  rows  (ii  ,i2)  and  (13,14)  can  be  connected.  If  there  is  a  pair  of 
treatments  (k,k)  appearing  in  any  row  (or,  column)  of  (j,j')  (or, 
(i,i'))  than  j~j'  (or,  i~i'). 

3.  A  Simple  Graphical  Method 

Suppose  the  pairs  of  treatments  (kj,k2)  az.d  (k3,k4)  are  occuring 
in  a  block  ( J 1 , J 2 )  and  the  rows  ij  and  i2  of  a  3  dimensional  design. 

It  follows  from  (1)  that 

E(yilJlkl  "  yilj2k2  "  yi2J  1^3  +  yi2j2^4)  “  TUi,k2)  "  T(k3,k4)-  (5 
We  therefore  say  that  (kj,k2)  and  (k3,k4)  are  connected  in  the 

analytical  sense  of  Definition  1  if  they  occur  in  a  block  (jj,j2),  or, 

in  other  words  (kj,k2)  and  (k3,k4)  are  connected  by  the  block  (jj,j2) 

which  contains  them.  From  Lemmas  1  and  2,  it  follows  that  (k^.k^)  and 

(k3,k4)  with  kj  -  k3 ,  occuring  in  two  different  blocks  may  be 

connected  because  some  other  pairs  of  treatments,  e.g.  (k5,k2)  and 


(k5,k4),  are  connected  by  some  other  block  containing  them.  In  this 
case,  we  say  that  (k3,k2)  and  (k3,k4)  are  connected  by  a  third  block. 

Suppose  (kj,k2)  is  occuring  in  a  block  C  j  i ,  j  2 )  an<*  a  row  *1  ar,d 
(k3,k4)  is  occuring  in  a  block  (j3,j4)  and  a  row  i2  of  a  3  dimensional 
design.  It  follows  from  (1)  that 


E(yiljlM-yil^2-yi2j3^3  +  yi2W 


(6) 


"  T(kltk2)  T(k3,kA)  +  6(jlfj2)  e(.13,j4)‘ 

If  (kj  ,k2)  and  (k3,k4)  are  connected  by  a  block,  then  it  follows  from 
(6)  that  C j 1 *j 2 )  and  ( j  3  » j  4 )  are  also  connected  and  vice  versa.  We 
say  combining  (5)  and  (6)  that  two  blocks  and  (j3,j4)  are 

connected  if  every  treatment  (kt  ,k2)  in  ( j 1 » j 2 )  is  connected  to  every 
treatment  (k3,k4)  in  ( j  3  , j 4 ) . 

We  now  introduce  a  new  concept  of  connected  path  for  three 
dimensional  designs.  Consider  a  graph  whose  points  are  the  pairs  of 
treatments  and  any  two  points  are  joined  by  an  edge  if  they  are 
connected,  in  the  sense  of  Definition  1,  by  a  block.  When  two  points 
(kj,k2)  and  (k3  ,k4 )  are  joined  by  a  block  (ji,j2)»  we  write 


(kj,k2)  (-W  (k3,k4). 


The  adjacency  of  any  two  points  in  the  graph  can  be  determined  from 
the  fact  that  all  points  occuring  in  the  same  block  of  the  design  are 
adjacent  and  applying  the  Lemmas  1  and  2. 


Definition  2.  Two  points  (kj,k2)  and  (k{  ,k'2 )  in  the  graph  are  said  to 
be  connected  if  there  is  a  path 

(ki,k2)  (jjiL2>  (k3,k4) - (kl,k2), 

joining  (kj,k2)  and  (k\,k2). 

It  follows  from  Theorem  1  that  the  design  is  connected  w.r.t. 
treatment  if  any  two  points  of  the  graph  are  connected  (i.e.,  the 
graph  is  connected).  Thus  for  a  design  connected  w.r.t.  treatment 
there  is  a  single  component  in  the  graph  (or,  equivalently,  there  is 
one  and  only  one  equivalence  in  K  as  discussed  in  Section  2). 

We  now  discuss  a  simple  procedure  of  determining  whether  a  design 
is  connected  or  not.  We  assume  without  any  loss  of  generality  that 
b  _<  a  and  consider  the  column  blocks  (1,2),  (1 ,3) ,...,( 1 ,b)  and 
(2,3).  Note  that  the  points  in  any  block  form  a  completely  connected 
subgraph.  Two  subgraphs  are  connected  when  any  point  of  one  subgraph 
can  be  joined  by  a  connected  path  to  any  point  of  the  other  subgraph. 
Lemmas  1  and  2  are  the  only  tools  needed  in  verifying  connectedness  or 
disconnectedness  of  two  subgraphs.  The  main  step  in  our  procedure  is 
to  check  whether  b  subgraphs  are  part  of  a  single  component  of  the 
graph  or,  in  other  words,  whether  these  b  subgraphs  are  connected  or 
not.  If  these  b  subgraphs  belong  to  different  components  of  the  graph 
(i.e.,  any  of  them  is  disconnected  from  some  others),  then  the  design 
is  not  connected  w.r.t.  treatment  and  block  (column) .  If  b  subgraphs 
belong  to  a  single  component  of  the  graph,  then  the  design  is 
connected  w.r.t.  block  (column).  This  in  fact  means  the  connectedness 
of  b(b-l)  subgraphs  corresponding  to  b(b-l)  column  blocks.  The 


V 


connectedness  of  the  subgraphs  (jj,12)  and  ( j  2  * J 1 )  implies  the 
connectedness  of  treatments  (k1,k2)  in  ( j  1 » j 2 )  and  (k2,kl)  in  (j2,j1) 
occuring  in  the  same  row.  This,  in  turn,  implies  the  connectedness  of 
two  treatments  kj  and  k2 ,  or,  in  other  words,  all  treatments  in  every 


row  of  the  design  are  connected. 


4 .  Examples 

We  now  present  two  examples  to  illustrate  the  procedure  discussed 
in  the  earlier  Section.  Our  first  example  of  a  3  dimensional  design 
is  taken  from  Shah  and  Khatri  (1973). 


1 

2 

5 

6 

3 

A 

7 

8 

8 

6 

1 

3 

7 

5 

2 

A 

The  column  blocks  (1,2),  (1,3),  (1,4)  and  (2,3)  are  shown  below. 


(1,2) 
1  2 
3  A 
8  6 
7  5 


(1.3) 
1  5 
3  7 
8  1 
7  2 


(1  ,9) 
1  6 
3  8 
8  3 
7  A 


(2,3) 
2  5 
A  7 
6  1 
5  2 


From  the  block  (1,4)  we  have 

(1,4)  (1,4) 

(1,6)  (3,8)  •  — -  (8,3). 

It  follows  from  Lemma  2  that  1~6  and  therefore  (8,6)  —  (8,1). 
Note  that  (8,6)  is  in  the  block  (1,2)  and  (8,1)  is  in  the  block 


,r 


(1.3) .  Thus  the  third  block  (1,4)  is  connecting  the  treatments  (8,6) 
and  (8,1).  This  fact  is  expressed  by 

0,4) 

(8,6)  (8,1). 

Similarly  from  the  block  (2,3), 

(2,3) 

(2.5)  (5,2), 

and  this  implies  2-5.  Thus  (1,2)  in  the  block  (1,2)  and  (1,5)  in 

the  block  (1,3)  is  connected  by  a  third  block  (2,3)  and  we  have 

(2.3) 

(1,2)  - - -  (1,5). 

It  is  interesting  to  observe 

(1.4) 
or 

,  x  (2,3) 

(1.6)  - - (6,1). 

The  subgraphs  corresponding  to  the  blocks  (1,2)  and  (1,3)  form  one 
connected  component  and  the  subgraphs  corresponding  to  the  blocks 

(1.4)  and  (2,3)  form  another  connected  component  in  the  graph. 
Therefore  the  design  is  not  connected  w.r.t.  treatment  and  block 
(column) . 


(7,5)  (8,6)  (8,1)  (7,2)  (8,3)  (7,4)  (4,7)  (5,2) 


- (2,3) - 

Figure  1.  Connected  and  disconnected  subgraphs  corresponding  to  the 
blocks  (1,2),  (1,3),  (1,4)  and  (2,3). 


The  second  example  of  a  3  dimensional  design  is  given  below. 


1 

2 

3 

2 

3 

4 

1 

3 

4 

1 

2 

1 

3 

(4,1). 

4  and  the  design  is 


From  the  column  blocks  (1,2)  and  (2,3)  we  have 

(1,2) 

or 

(1,2)  (2,3)  (3,4) 

From  Lemma  2,  it  follows  that  the  treatments  123 
completely  connected. 


5 .  Interrelation 

The  graphical  method,  BDS  algorithm  and  A  operator  in  Srivastava 

and  Anderson  are  techniques  in  solving  the  same  problem  of  determining 

connectedness.  There  is  a  common  ground  to  all  these  methods  and  this 

can  be  seen  from  the  equations  (5)  and  (6).  The  R-process  and  the 

Q-process  in  the  BDS  algorithm  can  be  applied  to  find  estimable 

constrasts  involving  only  r  parameters.  A  loop  L  defined  in  the  BDS 

algorithm  can  be  obtained  from  the  observation  y  ,  ,  y  ,  , 

ililkl  1 1-1 2^2 

y,  .  ,  and  y  ,  and  the  corresponding  T-functional  is  f(L)  = 

*2j2k4  i 2  J lk3 

T  ,  -  x  ,  .  ,  which  appears  in  (5).  Thus  the  observations 

(kl»k2)  (k3,k4) 

corresponding  to  the  treatments  (k,,k„)  and  (k„,k.)  in  the  rows  i,  and 
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±2  of  the  block  C j i » J 2 ^  gives  a  loop.  A  loop  can  also  be  a 
combination  of  several  such  loops.  A  quasi-loop  defined  in  the  BDS 
algorithm  is  in  fact  obtained  this  way.  From  the  connected  paths 
discussed  in  Section  3,  one  can  write  down  the  loops  and  the 
quasi-loops  in  the  BDS  algorithm.  There  is  obviously  an  even  number 
of  pairs,  in  the  sense  of  BDS  algorithm,  in  a  loop  or  a  quasi-loop. 
The  graphical  method  discussed  in  this  paper  is  lot  simpler  to  imple¬ 
ment  than  the  BDS  algorithm.  To  illustrate  this  let  us  consider 

Example  1  on  page  103  in  Birkes  et.  al.  (1976).  Assuming  y  *  t  and 

(2  3) 

considering  the  blocks  (2,3)  and  the  rows  2  and  3,  we  get  (2, A) 

(2  3) 

(2.1)  and  (4,4)  (2,1).  It  follows  from  Lemma  2  that  1-2  and 

1-4,  i.e.,  Tj  -  t2  and  t1  -  are  estimable.  Considering  the  block 

(12) 

(1.2)  and  the  rows  2  and  4,  we  find  the  treatments  (3,2)  (1,3). 

(2  3) 

Now  1-2  implies  (1,3)  *  (2,3).  Thus  we  have  the  connected  path 

(3.2)  ^  (1,3)  ^»3)  (2,3)  and  this,  in  turn,  implies  2-3  i.e., 

t2  “  t3  is  estimable.  We  do  not  need  to  write  the  matrices  N  and  M 
and  to  do  complicated  operations  as  in  the  BDS  algorithm. 

The  idea  of  a  loop  or  a  quasi-loop  is  present  in  Definiton  4.2 
and  in  the  condition  (a)  of  Theorem  4.2  in  Srivastava  and  Anderson 
(1970)  when  T-functional  is  a  simple  contrast  of  t's.  To  illutrate 
this,  suppose  (kj,k2)  and  (kj.kj)  are  occurring  in  (jj,j2)  and  the 
.  Then  T^  ■  (_to  ■  (i  1  ,j  1  ,ki  )  ,t?  m  (i2*j2>^l^)» 


rows  ij  and  i2 


T2  -  (li  -  (ii,j2.k2>.  13  -  d2  »j  1  »k2 ) ) »  6  -  2,  and  T  - 
(t^O »  ,  ^2 »  £3)  is  a  chain  connecting  the  levels  k^  and  k2  of 

treatment  in  the  notation  of  Srivastava  and  Anderson.  (Incidentally, 
there  are  some  misprints  in  the  equation  (4.6)  and  the  other  places, 
although  the  results  are  all  correct.) 


6.  Higher  Dimensional  Designs 

The  method  discussed  in  Sections  2-5  for  3  dimensional  designs 
can  be  generalized  for  designs  of  dimensions  4  or  more.  Consider  a  4 
dimensional  design  with  a  blocks  of  type  I,  b  blocks  of  type  II,  c 
blocks  of  type  III  and  v  treatments.  The  standard  additive  model  is 
E*yijJlk>  “  ^  +  ■*"  "*■  T  ^ »  i“l,...,a,  j*l »...  ,b, 

£*l,...,c,  k=l , . . . ,v.  (7) 

Suppose  the  treatments  kj ,k2 ,k3  and  k4  are  occurring  in  the  blocks  i j 
and  i2  of  type  I,  ^  and  j2  of  type  II  and  £j ,  £2 ,  £3,  £4  of  type 
III.  Then 


E(y 


il J lAlkl  yil32£2k2  ri2Jl£3k3 


-  y. 


^1 2 J  2^4k4  ^ 


( (k4  ,k3 )  ,  (k2,kj))  ( ( ^4 » *3 ) ,  (^2»^i))» 

whe  re  x  s  x  ~  t 

((k4,k3),  (k2  ,k^ ))  (k2,k})  (k4,k3), 

’  Y(*2.*i)  "  Y(*4,*3)  and  Y(*2,£i)  *  Y«i  "  YV 

design  is  thus  reduced  to  a  2  dimensional  design. 


(8) 

Y((*4,*3>.  (*2»*l)> 

A  4  dimensional 

Any  two  treatments 


occurring  in  a  block,  say  ((£4,£3),  (£2,£i)),  are  connected.  Any  two 
blocks  are  connected  if  treatment,  ((ki+,k3),  (k2,kl)),  in  one  block, 


-14- 


say  ((£4,^3),  (&2»*l))»  is  connected  to  a  treatment,  say  ((k^.k^), 
(k^.kj)),  in  the  other  block,  say  ((£4,^3),  The  connected¬ 

ness  or  disconectedness  w.r.t.  treatment  and  type  III  block  can  easily 
be  determined  by  defining  an  equivalence  relation  on  appropriate 
sets.  The  determination  of  connectedness  w.r.t.  type  I  block  and  type 
II  block  is  then  straightforward.  The  idea  of  reduction  in  dimension¬ 
ality  jutifies  the  partitioning  of  the  parameters  in  the  model  (7) 
using  the  BDS  algorithm,  see  Example  2  of  Section  10  in  Birkes,  Dodge, 
and  Seely  (1976). 


Remarks 

If  the  model  is  nonadditive,  i.e.,  interactions  are  present  in 
the  model,  then  the  idea  of  reduction  of  dimensionality  does  not 
work.  However,  the  idea  in  this  paper  helps  in  identifying  the 
estimable  parametric  functions.  To  illustrate  this,  we  consider  a 
nonadditive  model  for  a  3  dimensional  design 


E(yyk>  -  u  +  +  UB)^  +  V  (8) 

where  (aB)^  is  the  interaction  effect  between  the  I1"*1  row  and  the  jth 


column.  A  loop  can  still  be  obtained  from  the  observations  y  , 

1 1 J  lk  1 

yj  .  ,  ,  y,  .  .  and  y,  ,  and  the  corresponding  t,  (aB)  - 

*-2  j  2*4  a2J  lk3 

functional  is  x,  -  x,  -  x,  +  x,  +  (aB),  ,  -  (a&),  ,  - 

*1  k2  k3  k4  i l j l  i l J  2 

(aB),  ,  +  (aB),  ,  ,  the  right  hand  side  of  (5)  under  (8). 

i2jl  12  j  2 
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